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On the existence of chaotic circumferential waves in spinning disks 
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We use a third-order perturbation theory and Melnikov's method to prove the existence of chaos 
in spinning circular disks subject to a lateral point load. We show that the emergence of transverse 
homoclinic and heteroclinic points respectively lead to a random reversal in the traveling direction 
of circumferential waves and a random phase shift of magnitude 7r for both forward and backward 
wave components. These long-term phenomena occur in imperfect low-speed disks sufficiently far 
from fundamental resonances. 
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Transversal vibration modes of hard disk drives 
(HDDs) are excited by the lateral aerodynamic 
force of the magnetic head. Previous works (J, [2[ 
revealed that chaotic orbits are inevitable ingre- 
dients of phase space flows when the lateral force 
is large, or the disk is rotating near the critical 
resonant speed. For low-speed disks, however, an 
adiabatic invariant (a first integral) was found [2] 
using a first-order averaging based on canonical 
Lie transforms. According to the first-order the- 
ory, regular vibrating modes of imperfect, low- 
speed disks are independent of the angular veloc- 
ity of the disk, Qq. In such a circumstance, the 
speed of circumferential waves is the natural fre- 
quency of the lateral mode, u>, derived from linear 
vibration analysis. HDDs are usually operated 
with angular velocities smaller than lu (safely be- 
low resonance). Moreover, the magnitude of the 
lateral force F is very small. Given these condi- 
tions, we show that it is impossible to continue 
the Lie perturbation scheme [2] up to terms of 
arbitrary order and remove the time variable t 
from the Hamiltonian. In fact, due to the spe- 
cial forms of nonlinearities in the dynamical equa- 
tions of spinning disks, one can not remove t from 
third-order terms. Subsequent application of a 
second-order Melnikov theory reveals that trans- 
verse homoclinic and heteroclinic points do exist 
for all F, Qq ^ 0. This implies chaos, or equiva- 
lently, non-integrability of governing equations. 



I. INTRODUCTION 

Dynamics of continuum media, like fluids, rods, plates 
and shells is usually formulated as a system of partial 
differential equations (PDEs) for physical quantities in 
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terms of the spatial coordinates x and the time t as 



C(u) 



0. 



(1) 



Here £ is a nonlinear operator and u(x, t) is the vector 
of dependent variables. When the boundary conditions 
are somehow simple, approximate variational methods 
based on modal decomposition and Galerkin's projection 
[H can be used to reduce the order of governing equa- 
tions. These methods begin with solving an auxiliary 
eigenvalue problem, which is usually the variational field 
equation <$£(u, A) = 0, and build some complete basis 
set Uk(x, Ak) for expanding u(x, t) in the spatial do- 
main. Here Ak is an eigenvalue, which is characterized 
by the vectorial index k. Each Uk is called an eigenmode 
or a shape function. Such a basis set should preferably 
satisfy boundary conditions and be orthogonal. 

Once a complete basis set is constructed, one may sup- 
pose a solution of the form 



u(x,<) = ^Vk(t)-U k (x,A k ) 



(2) 



Substituting from © into (P) and taking the inner prod- 
uct 



J £(u)-U k '(x)dx = 0, 



(3) 



leaves us with a system of nonlinear ordinary differen- 
tial equations (ODEs) for the amplitude functions V"k(t). 
The evolution of the reduced ODEs shows the interaction 
of different modes and their influence on the development 
of spatiotemporal patterns. 

In a series of papers, Raman and Mote [l], Q used the 
modal decomposition method to investigate transversal 
oscillations of spinning disks whose deformation field is 
described in terms of the displacement vector (u,v,w) 
with u and v being the in-plane components. The most 
important application of the spinning disk problem is in 
design, fabrication and control of HDDs. The governing 
PDEs for the evolution of displacement components were 
first derived by Nowinski [f| and reformulated more re- 
cently by Baddour and Zu @. Let us define {R,qb) as 
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the usual polar coordinates. For a rotating disk with the 
angular velocity il , Nowinski's theory assumes that the 
in-plane inertias (u,u)Oq, 2Q,o(u t t,v t t) and {u it t,v,tt) are 
ignorable against RQq. This is a rough approximation for 
high-speed disks and one needs to use the complete set 
of equations as Baddour and Zu [6| suggest. Nowinski's 
theory, however, has its own advantages (like the exis- 
tence of a stress function) that facilitate the study of the 
most important transversal modes. In low-speed disks, or 
disks with high fiexural rigidity, one has il <C uj. Hence, 
it is legitimate to apply Nowinski's governing equations 
in such systems. 

In this paper we analytically prove the existence of 
chaos, and therefore, non-integrability of the reduced 
ODEs that govern the double-mode oscillations of im- 
perfect spinning disks. We investigate low-speed disks 
subject to a lateral point force exerted by the magnetic 
head. The lateral force in HDDs is very small and its 
origin is the aerodynamic force due to air flow in the gap 
between the disk and the head. We show that chaotic 
circumferential waves dominate some zones of the phase 
space over the time scale t ~ 0(e~ 3 ) with e being a small 
perturbation parameter. This indicates very slow evolu- 
tion of random patterns, and the practical difficulties of 
their identification. 

The paper is organized as follows. In ETTA we present 
the Hamiltonian function in terms of Deprit's [7| Lis- 
sajous variables. In £1III[ we use a canonical perturba- 
tion theory @, Q to eliminate the fast anomaly I from 
the Hamiltonian. The action associated with I then be- 
comes an adiabatic invariant. Transversal intersections 
of destroyed invariant manifolds, and therefore, non- 
integrability of the normalized equations, is proved by 
a second-order Melnikov method in ^IVl We present a 
complete classification of circumferential waves in EfVland 
end up the paper with concluding remarks in 3 VII 



II. PROBLEM FORMULATION 

Let us assume U m (R) as an orthogonal basis set that 
represents the disk deformation in the radial direction. 
The index m stands for the number of radial nodes 
that U m (R) has. According to Raman and Mote's [l| 
treatment of imperfect disks, the following choice of the 
transversal displacement field 

w(R,4>,t) — U m (R) [x (t) cos ncf) + y (t) sin n<f>], (4) 

reduces Nowinski's governing equations to a system of 
ODEs for the amplitude functions x(t) and y(t) as 

x + \ 2 x + e-f (x 2 + y 2 ) x = eF cos(nf2 £) , (5a) 
y + uj 2 y + £7 (x 2 + y 2 ) y — eFsm(nflot), (5b) 

where A, to and 7 are constant parameters that depend 
on the geometry and material of the disk, e is a small 
perturbation parameter, F is the weighted integral of the 



lateral point force, and f^o is the angular velocity of the 
disk. 

We suppose small deviations from perfect disks and 
write the constant parameter of ([5]) as X 2 /uj 2 = 1 + e-q. 
We also define nQ,$ — eO with 0(Q) ~ 0(u>). Denoting 
(jPxiPy) as the momenta conjugate to (x,y), it can be 
verified that equations ([5]) are derivable from the Hamil- 
tonian function 



H 



(pI+P 2 v) + ^ 2 (x 2 +y 2 )+> 



np+- (x 2 

4 v 



1 



V 2 2 
—CO X 

2 



y 2 ) 2 — F (xcosp 



ysinp) . (6) 



We have introduced the action P and its conjugate an- 
gle p = eQt to make our equations autonomous, which 
is a preferred form for the application of canonical per- 
turbation theories. The extended phase space has now 
dimension six. Dynamics generated by ([6]) is better un- 
derstood after carrying out a canonical transformation 
(x,y,p x ,p y ) — > (l,g,L,G) to the space of Lissajous vari- 
ables Q so that 



x = scos(g + I) — dco$(g — I), 

y = s sin(g + I) — dsin((? — I), 

p x = — uj [s sin(<7 + I) + dsm(g — I)] . 

p y = uj [scos(g + I) + dcos(g — I)] , 



(7a) 
(7b) 
(7c) 
(7d) 



G 



2uj 



L-G 

2uo ' 



L>0, \G\<L. 



In the space of Lissajous variables, the Hamiltonian de- 
fined in ^ becomes 



H 


= H (L) + 


eH.il, g,p,L,G,P), 


(8) 




= luL, 






Hi 


= np-F 


s cos(g + I — p) — G?cos(g 
d 2 ) - 2sdcos(2Z)] 2 


-l-p)} 




+ ^ 


+ d 2 ) + s 2 cos(2.g + 21) - 


2sd cos(2i?) 



d 2 cos(2g - 21) - 2sdcos(2l) 



From ([8]) we conclude that I is the fast angle, and g and 
p are the slow ones. Therefore the long-term behavior 
of the flows generated by ((5|) can be analyzed by averag- 
ing H over I. After removing I, its corresponding action 
L will be a constant of motion for the flows generated 
by the averaged Hamiltonian (H)i, and the phase space 
dimension reduces from 6 to 4. 



III. CANONICAL THIRD-ORDER AVERAGING 

In order to average H over I, we use the normal- 
ization procedure of Deprit and Elipe ^9j. Denoting 
X = (l,g,p) and Y = (L, G, P), we define a Lie transfor- 
mation (l,g,p,L,G,P) — > (l,g,p,L,G,P) as 



X = E W (X), Y = E W {Y), 



(9) 
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so that the Hamiltonian function in terms of the new 
variables, K = (H)t, does not depend on I. E w is the Lie 
transform generated by the function W and it is defined 

as 



E W (Z) = Z+(Z;W) + ^((Z;W);W) 
+ l(((Z;Wy,W);W) + --- . 



(10) 



In this equation, (/i;/2) denotes the Poisson bracket of 
/i and f 2 over the (l, g,p, L, G, P)-space. We expand the 
generating function W as 



W = eWi 



2! 



W 2 



(11) 



and specify the averaged, target Hamiltonian K = 
K(g,p, L, G, P) as the series [1] 



K = K + eiiTi 



2! 



K 2 



3! 



K 3 + 



with 



K = uiL, 
1 



oZ / M ' 
2tt Jo 



(12) 

(13a) 
(13b) 



K 2 = f- / [2{H X] Wi) + ((H ; W x )\ Wx)]dl, (13c) 
K 3 
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/ \3(H 1 ;W2)+3((H 1 ;W 1 );W 1 ) 

2((# ; W 2 ); Wk) + ((H ; Wi); VK 2 ) 



(((#o;Wi);Wi);Wi) 



dl. 



(13d) 



Wi and W2 are determined through solving the following 
differential equations 



dW 1 

~df 

dW 2 

W 



H^l, g,p, L, G, P) — Ki(g,p, Z, G, P),(14a) 
2(H l -W 1 ) + ((H -Wi)-W l )-K 2 . (14b) 



By substituting from (TTJ)) into (fl3j) and evaluating the 
integrals, one finds the explicit form of the new Hamilto- 
nian K, which has been given in Appendix [A] up to the 
third-order terms. 

Once I is removed from the Hamiltonian, L becomes 
an integral of motion. The slow dynamics of the system 
is thus governed by the flows in the (g, G)-space. We 
introduce the slow time r = p/fl, ignore the fourth-order 
terms in e, and obtain the following differential equations 
for the dynamics of (<?, G) 



dg _ dK 
dZ~dG~ 
dG _ dK 
dr dg 



fiti,G)+ehifaG,e,T), (15a) 
= f 2 (g,G) + eh 2 (g,G,e,T), (15b) 



where 

fi = d 6 + di cos(2g), f 2 = ei sin(2g), 
hi = dj + d 2 cos(2g) 

+ e[d$ + d 3 cos(2g) + d± cos(4g) + d§ cos(2g — 2f2r)], 
h 2 = e 2 sin(2g) 

+ e[e 3 sin(2<?) + e 4 sin(4c7) + e 5 sin(2<? - 2Q,t)}. (16) 



In these equations, di [i = 1, • 



ande, (j = !,••• ,5) 



are functions of Z and G (see Appendix |A"|) . It is re- 
marked that the action P appears only in K\ via the 
term flP. It then disappears in the normalized equa- 
tions (fTS")) after taking the partial derivatives of K with 
respect to g and G. The partial derivative of K with 
respect to P determines the evolution of p, which is in 
accordance with the simple linear law p(r) = Qt + p(0). 
The dynamics of P itself is governed by 



d_P 
"d7 



dK 
dp 



l_d_ 

ndZ 



K(g,G,r) 



(17) 



One may integrate (fl~7|) to obtain P(t) once equations 
(fTof are solved. The behavior of P is thus inherited from 
g(r) and G(r). 



IV. THE MELNIKOV FUNCTION 

There are few analytical methods in the literature for 
the detection of chaos in perturbed Hamiltonian systems 
[ToL [TlT |. Melnikov's [l(| method is the most powerful 
technique when the governing equations take the form 



dx 
d7 



f(x) + eh(x, e, r), x G 



(18) 



so that the unperturbed system dx/dr = f(x) is in- 
tegrable and possesses a homoclinic (heteroclinic) orbit 
q h (r) to a hyperbolic saddle point, and h(x, e,r) is T- 
periodic in r. The occurrence of chaos is examined by 
the Melnikov function 



M(r , e) = eMi(7ft) + e 2 M 2 (r ) + 



(19) 



where Mfc(ro) denotes the fcth-order Melnikov function. 
Assume that Mi(ro) is the first nonzero term, i.e., 
M fe (r ) = for 1 < k < i - 1. If M^tq) has simple 
zeros, then, for sufficiently small e, the system (|T8|) has 
transverse homoclinic (heteroclinic) orbits, which i mply 
chaos due to the Smale-Birkhoff homoclinic theorem [12| • 
The first-order term in (fH?)) is determined by the classical 
formula 



Mi (t ) 



-Hoc 



f (q h (t)) A h (q h (r) , 0, r + r ) dr, (20) 



where the wedge operator A is defined as f A h 
fih 2 — j 2 h\. Although the Hamiltonian equations fTE 
have a suitable form for the application of Melnikov's 
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method, they are autonomous up to the first-order terms 
in e. Consequently, M 1 (r ) vanishes identically for all 
t G [0, T}. We thus need to investigate the second-order 
Mclnikov function. For doing so, we begin with solving 
the unperturbed system 



d| 
dr 
dG 
dr" 



dKi 
~dG ~ 
dK 1 
dg 



= / 2 (.g,G), 



(21a) 
(21b) 



along homoclinic (heteroclinic) orbits. Jalali and An- 
goshtari showed that for L > t\lo 3 j^f, equations (12"T]) 
have hyperbolic stationary points at Sq = (go, Go) — 
(-tt.O), S 2 = (92, G 2 ) = (0,0), and 5 4 = (g 4 , G 4 ) = 
(■n, 0). The implicit equation of the invariant manifolds 
that terminate at the saddle points are 



cos[2g h (i 



L 



■Gl(r) 



[L 2 Gl(r)] 



-1/2 



(22) 



For 7Z > 2t]u> 3 , equation (|22|) represents a heteroclinic 
orbit which connects 6*0 to S 2 - For ijuj 3 < 7Z < 2-qui 3 , 
the heteroclinic orbit disappears and it is replaced by a 
homoclinic orbit (see Figure [l}. To compute the explicit 
form of the homoclinic (or heteroclinic) orbit of ([2Tj) , we 
use (|21b|) and ([22]), and obtain 



Gh(r) 



dG 



G h (0) GV 1 - aG 2 

J 7Z — f?CJ 3 



(23) 



7 

16?y/3w 4 



4w 



where the lower integration limit is Gh(0) = l/\/a. After 
taking the integral (|23p . we arrive at 



G h (r) = 



sech(v73 r) 



(24) 



for the G > branch of the homoclinic (heteroclinic) 
orbit. Having G^(t), it is straightforward to calculate 
cos[2gh(r)] and sinp^hML an d determine the explicit 
form of q h (r) = [g h (r), G h (r)] . 

For constructing Af 2 (ro), we use Franchise's [1J, [l4| 
algorithm that has been devised for dynamical systems 
with polynomial nonlinearities. To express the averaged 
Hamiltonian K in terms of polynomial functions of some 
new dependent variables, we utilize Hopf's variables 



Ql = -L^ 2 -G 2 cos(2g), 



2w 



0.2 = ^VlI-G* sin(25), 
2ui 



(25a) 
(25b) 



and obtain the following differential 1-form for the evo- 
lution of the averaged system 

0Q1 oQ 2 
+ e [(si + es 2 ) dQ 2 - ( Zl + ez 2 ) dQi] = 0. (26) 



Here, the first-order Hamiltonian is 



Ki(Q 1 ,Q 2 ) = ^(Q 2 1 + Ql)-^-Qi 

and 



G. 



(27) 



si = miQ 2 , (28a) 
zi = riiQi + n 6 , (28b) 
s 2 = m 2 (Q\ + Q 2 )Q 2 + m 3 QiQ 2 + m, i Q 2 

+ m 5 sin(2f>r), (28c) 
z 2 = n 2 (Q\ + Q 2 2 )Q 1 + n^Ql + Ql) + n 4 Qi 

+ n 5 cos(2ilr) + n 7 . (28d) 

The constant coefficients G, m.; (i = 1, • • • ,5), and rij 
(j = 1, • • • ,7) have been given in Appendix IBI A prereq- 
uisite for the application of Franchise's [l3| algorithm is 
that for all polynomial 1-forms D that satisfy the condi- 
tion 



D = 0, 



(29) 



there must exist polynomials A(Qi,Q 2 ) and r(Qi,Q 2 ) 
such that D = dA + rdKi. We call this the condition (*) 
and prove in Appendix [Cl that K\ satisfies the condition 

(*)• 

Franchise's algorithm states that if M\(to) = ■■■ — 
A/fe_i(To) = for some integer k > 2, it follows that 



M fc (ro) = / D k , 
Jib 

D 1 = 8 X , D m 



E 

i-\-j—m 



nSj, 



5j = ZjdQi - SjdQ 2 , 



(30a) 
(30b) 
(30c) 



for 2 < m < k. The functions 77 are then determined 
successively from the formulas Di = dAi + TidK\ for 
i = 1, ■ ■ ■ , k — 1. We have already found that 

Mi (r ) = f <S 1 = 0, (31) 

Si = (m<2i + n 6 )dQi - m,iQ 2 dQ 2 . (32) 
From (f3"0")) and (IC2[) it can be shown that 



M 2 (t ) 



Substituting (|3"0c]) and (gSJ) into 
the integration along q^i"), result in 



(33) 

, and carrying out 



(34) 



M 2 (to) = ^-^sin(2r!ro), 

with / being a constant (see Appendix [D|) . Equation 
shows that r„ = nir/ (2il) (n = 1, 2, • ■ • ) are simple zeros 
of M 2 (to) so that 



M a (r„) = 0, 



dM 2 (T Q ) 



8t 



(35) 
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FIG. 1: Possible topologies of the phase space flows of the averaged system for F = 0. In all panels the horizontal axis indicates 
g/n and the vertical axis indicates G/L. In the first topology (left panel), all stationary points are of center type. In the second 
topology (middle panel) two new centers, with non-zero G-coordinates, have emerged for g — ±nir (n = 0, 1) and symmetrical 
homoclinic loops (thick lines) connect saddle points to themselves. In the third topology (right panel) the off-axis centers (and 
their surrounding tori) are still present but the separatrix curves (thick lines) are of heteroclinic type. 




FIG. 2: Phase space structure of the normalized equations (|15[) with L = 2, u> — rj = F — 1, and Q, — e — 0.1. Left panel: 
7 = 0.8. Right panel: 7 = 2. In both panels the horizontal and vertical axes indicate g/n and G/L, respectively. 



Thus, we conclude that the global stable and unstable 
manifolds of the saddle point S£° , W S (S;<> ) and W U (S™), 
always intersect transversely. Transversal intersections 
cause a sensitive dependence on initial conditions due to 
the Smale-Birkhoff homoclinic theorem. This is a route 
to chaos. On the other hand this means that the reduced 
equations (fT5)) are non-integrable for F 7^ 0. 



V. CLASSIFICATION OF CIRCUMFERENTIAL 
WAVES 

For F = 0, h does not depend on t and the normalized 
equations (|15p are integrable. In such a circumstance, the 
phase space structure can take three general topologies 
(depending on the values of the system parameters and 
L) as shown in Figure [TJ In the first topology all sta- 



tionary points with the coordinates (<? S ,G S ) calculated 
from 

f(g s ,G s ) + eh(g s ,G s ,e) = 0, (36) 

are centers and they lie on the G = axis with g s = 
—ir + mr/2 (n = 0, • • • ,4). In the second and third 
topologies, two off-axis centers (with G s 7^ 0) come to 
existence for g s = ±nn [n = 0,1) and the on- axis sta- 
tionary points with the same g s = ±nir become saddles. 
In the second topology, each saddle point is connected to 
itself by a homoclinic orbit, and in the third topology, a 
heteroclinic orbit connects two neighboring saddle points. 
The system with heteroclinic orbits allows for rotational 
g(r) while in the system with homoclinic orbits g(r) is 
always librating. Beware that this classification of phase 
space flows is valid as long as e is sufficiently small. 
For F = 0, the phase space flows of (p~5|) are 
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structurally stable (with no unbounded branches) and 
the whole (g, G)-space is occupied by periodic orbits 
of period T(K). At the stationary points, one has 
T(K(g s ,G s ,L)) — 0. Given the invariance of L, and 
the periodic solutions g(r) = g(r + T(K)) and G(r) = 
G(t + T(K)), the anomaly I is determined through solv- 



ing 



— = lu + e— - K x + —K 2 + —;K 3 
dt dL V 2! 3! 



(37) 



which results in I = cut + eR(r) with R(t) = R(t + 
T{K)). According to ©, the functions g(t), G(t) and 
L(t) are also periodic in t and we conclude that l(t) = 
ut + Rw(t) with R\v(t) = Ew(l) — being a small- 
amplitude periodic function of t. The explicit from of 
the circumferential wave will then become 



w(R,</>,t) _ L(t) + G (t) 
U m (R) V 2to 



■ cos [n<f>— cot— Rw{t)— g(t)] 



m-G(t) 



2uj 



cos [ncf)+u)t+R w {t)-g{t)} , (38) 



which is composed of a forward and a backward traveling 
wave. Due to the periodic nature of L{t) and G(t), when 
the amplitude of the forward traveling wave is maximum, 
that of the backward wave is minimum and vice versa. 

As our results of ^1 VI shows, the regular nature of trav- 
eling waves is destroyed for F ^ and a chaotic layer oc- 
curs through the destruction of the homoclinic and hete- 
roclinic orbits of (pH]) . This happens over the time scale 
t ~ 0(e~ 2 ) or t ~ 0(e~ 3 ) (because p is present only 
in K3). Figure [2] shows Poincare maps of the system 
(|15p for F =/= 0. The sampling time step in generating 
the Poincare maps has been 2ir/tt. It is seen that most 
tori around elliptic fixed points are preserved. They cor- 
respond to regular periodic and quasi-periodic solutions 
of the normalized system. For chaotic flows, the func- 
tions g(r) and G(r) randomly change within the invari- 
ant measure of the chaotic set. Consequently, the orig- 
inal Lissajous variables g(t), L(t), G(t) and also Rwif) 
become chaotic too. 

For G(t) > and G(r) < the forward and the back- 
ward traveling waves are the dominant components of 
the circumferential wave, respectively. When the chaotic 
layer emerges from the destroyed homoclinic orbits (left 
panel in Figured]), the sign of G(r) is randomly switched 
along a chaotic trajectory. This means a random trans- 
fer of kinetic/potential energy between the forward and 
backward traveling wave components. For chaotic trajec- 
tories of this kind the angle g(r) randomly fluctuates near 
g w ±mt (n = 0, 1) with an almost zero average. The 
evolution of circumferential waves is quite different when 
the chaotic layer emerges due to destroyed heteroclinic 
orbits (right panel in Figured]). In this case chaos means 
a random change between the librational and rotational 
states of g(r). Such a change induces an unpredictable 
phase shift of magnitude tt for both forward and back- 
ward traveling wave components. We note that G(r) can 



flip sign on a chaotic trajectory only when g(r) is in its 
librational state. 



VI. CONCLUDING REMARKS 

Resonance overlapping [TJ [l6[ is the main cause for 
chaotic behavior in spinning disks with near-resonant an- 
gular velocities The chaos predicted in this paper, 
however, happens far from fundamental resonances. Op- 
tical and HHDs are usually operated below critical res- 
onant speeds and the lateral force F due to magnetic 
head is very small. We showed that whatever the magni- 
tude of F may be, a chaotic layer fills some parts of the 
phase space because the Melnikov function of the nor- 
malized equations has always simple zeros. Dynamics of 
rotating disks is regular only if F vanishes, which is an 
unrealistic assumption for disk drives. In low-speed disks 
with small F, diffusion of chaotic orbits (within their in- 
variant measure) takes a long time of t ~ 0(e~ 3 ). The 
slow development of chaotic circumferential waves makes 
them undetectable in short time scales at which most 
controllers work. The Melnikov function (|34[) depends 
not only on F, but also on the parameter rj through the 
constant I. The parameter 77 is a contribution of imper- 
fections, which are likely because of limited fabrication 
precision in micro/nano scales. For a perfect disk with 
rj = 0, the off-axis elliptic stationary points of (f!H]) . and 
consequently, homoclinic and heteroclinic orbits disap- 
pear. In such a condition the Melnikov function is indef- 
inite, but the system admits an exact first integral and 
the dynamics is governed by the Hamiltonian function 
given in equation (11) of Jalali and Angoshtari 0]. 

One of the most important achievements of this work 
was to unveil the fact that it is premature to truncate the 
series of canonical perturbation theories before recording 
the role of all participating variables. In systems with 
non-autonomous governing ODEs (non-conservative sys- 
tems), one must be cautious while removing a fast angle 
through averaging schemes. The removal of the fast an- 
gle may also wipe out time-dependent terms, up to some 
finite orders of e, and hide some essential information of 
the underlying dynamical process. Strange irregular so- 
lutions can indeed occur at any order and influence the 
long term response of dynamical systems as we observed 
for the spinning disk problem by keeping the third-order 
terms. 
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APPENDIX A: THE NORMALIZED 
HAMILTONIAN 



By evaluating the integrals in (|T3|) , we obtain the first, 
second and third order terms of the normalized Hamilto- 



Defining S = (L 2 - G 2 )/G, the functions ej{L,G) in 
equations (HHJ) become 

ej = -2djS, j = 1,...,5, j ^ 3, 



e3 = 



16r7 3 w 6 



128w 5 
9QL-fr] 2 ui 3 ] \fl? 



17(10L 2 - G 2 ) 7 2 ?7 



G 2 . 



(A5) 



-7G 2 L(3 7 L + 2r7w 3 ) 



OP 



Kn 



8u 2 8uj 2 
IJVZ 2 -G 2 cos(2g), 

-— ^ [2Z(-9G 2 + 17Z 2 ) 7 2 + 64F 2 w 3 
8(3Z 2 - G 2 ) 7 77w 3 + 8Zr7 2 w 6 
4w 3 [-6i7?/ - 2?/ 2 w 3 ] VZ 2 - G 2 cos(2g) 
3 



(Al) 



(A2) 



:{llG 4 7 3 - 258G 2 Z 2 7 3 + 375L 4 7 3 



512w 8 

1024F 2 Z 7 w 3 - 180G 2 Z7 2 t?w 3 + 340Z 3 7 2 r;w 3 
256F 2 nuj 6 - 48G 2 7?7 2 w 6 + 176Z 2 777 2 w 6 



32L7 ? 3 w 9 - 2uj- 



17(10i 2 - G 2 )7 2 ?7 



+ VQL-frfuj 6 + I677V 



VZ 2 - G 2 cos(2ff) 
16(G 2 - Z 2 ) 7 77 2 w 6 cos(45) + 256F 2 nuj 6 cos(2p) 
512F 2 7 w 3 VZ 2 - G 2 cos(2g - 2p)} . (A3) 



Consequently, the functions di(L,G) in equations (TT 
are found to be 



* = !^(Z 2 -G 2 )- 1/2 



G 
~16w 2 
3G 



(I4 
(h 
de 
d r 
da 



256w 5 
+ 96Z 7 ?7 2 u; 3 

3G 7 r7 2 



(6L 1V + 2r ] W)(L 2 -G 2 )-^ 2 , 

16?7 3 u; 6 + 51(4Z 2 -G 2 )7 2 ?7 
(Z 2 -G 2 )- 1 / 2 , 



_16w 2 
3G^ 2 7 

7G 
4w 2 ' 

G 



(Z 2 -G 2 )- 1 / 2 , 



32u}' 



(I8L7 2 + 87?7w 3 ), 



3G 



256w 8 



27 3 (11G^ - 129L' Z ) - 180Lj 2 r]Oj 3 



487?7 2 a; 6 



(A4) 



APPENDIX B 

The constant coefficients of equations (f27f and Ip8|) are 
as follows 



G = 

Til = 

"2 = 

ri3 = 

■124 — 

n 5 = 

n 6 = 

n 7 = 

m l = 

m 3 = 

7TI4 = 



L{~/L + t]lu 3 ) 
AJ 2 ' 

^(9Z 7 2 + 4777a; 3 ) 

33 7 3 
_ 48w 4 ' 
17 7 2 77 

1 



32w 6 
F 2 7 



(59L 2 7 3 + 45L7 2 ryw 3 + 1677? 2 w b ), 



2, ,3 



3L777 + ?7' u> 



8uj 



1 



(153Z 2 7 2 ?y + 96L 7 ?7^ 3 + 16^^) 



2, .3 



„3, ,6\ 



256w 4 

* — 1, 2, 5, 

-2n 3 , 

^-^(59Z 2 7 3 + 45Z 7 5 V 3 + 8 7 77 2 w 6 ). 



APPENDIX C 



(Bl) 



In this appendix we prove that Ki{Qi 1 Q 2 ) given in 
(f27|) . satisfies the condition (*). To this end, we need the 
following theorem. 

Theorem 1. Any polynomial 1-form D of degree n in 
Qi and Q 2 can be expressed as 

D = dA + rdKi + €(K 1 )Q 2 dQ 1 , (CI) 

where ^4(Qi, Q 2 ) and r(Qx, Q2) are polynomials of degree 
(n+1) and (n—1) respectively, and^(Ki) is a polynomial 
of degree [h(n— 1)] where [x] denotes the greatest integer 
in x. 

Iliev [l5[ has proved the same theorem for H = {Q\ + 
Theorem 1 can thus be proved in a similar man- 
ner. Here we only present a useful result. 

Let D be a general polynomial 1-form of degree 1, 



D = (oioQi + a iQ 2 + aoo)dQi 
+(b w Qi + b iQ 2 + b 00 )dQ 2 , 



(C2a) 
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then in (|C1[) we have 

KQi,Q 3 ) = 0, 

^(ATi) = a i - b 10 . 



(C2b) 
(C2c) 
(C2d) 



Since dlCi = along any phase space orbit characterized 
by Ki(Qi, Q2) — k, and since the integral of an exact 
differential dA around any closed curve is zero, from (|C1|) 
we obtain 

f D = e(A) / Q 2 dQi 
On the other hand, from (125aj) we have 



dQ 



dr 



i = £7(r) - 2Q 2 , 



where -E(r) is an even function of r. Given the fact that 



Q2 is an odd function of r, we conclude that 



/ D = -2£(fc) / Q2 d7 

Jq h J-oc 



Consequently, if J q D = 0, it follows that £(k) = and 
therefore Z? = dA + rdKi , which completes the proof. 



APPENDIX D 

In equation (|34| . the constant coefficient 7 is 

64tt^^ 2 w 2 , /vrfi\ 
I = 1 csch — = 

nUy/rj(3yL - 4?7w 3 )(/3 + 4ft 2 ) ( nfl 



7/3 

™LJr}(ZP - An 2 ) , /vrft 
1 sech — — 



sech I — — 
7 VV^ 
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